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THEORETICAL ANALYSIS OF SEVERAL LINEAR STEADY-ELOW 


PLASMA ACCELERATORS WITH CROSSED ELECTRIC 

AND MAGNETIC FIELDS 

By Clarence W. Matthews 
Langley Research Center 

SUMMARY 


The ideal one -dimensional equations which express the motion of an invis- 
cid, compressible, conducting, single-fluid plasma through crossed electric and 
magnetic fields are written in a nondimen sional form. Solutions of these equa- 
tions are then presented for conditions wherein the cross-sectional area is 
held constant or is allowed to vary in such a way that the velocity, pressure, 
or density is held constant. These conditions are considered for a number of 
variations of the current density with velocity. Several special cases are 
also presented in which the velocity, pressure, temperature, or Mach number is 
held constant in a constant-area channel. 

For illustrative purposes computed results are presented for most of the 
solutions. These results are used to show the nature of the various modes of 
channel operation and may be used to compare modes and to make a preliminary 
selection of the mode of operation required to achieve a given purpose. These 
computations show whether a channel will accelerate or decelerate the flow, 
whether the flow is regular or forbidden, whether the velocity is limited or 
unlimited as the channel length becomes infinite, or whether other unexpected 
peculiarities exist. 


INTRODUCTION 


The need for high-velocity gas flow has in the past several years promoted 
studies of the use of crossed electric and magnetic fields to accelerate a con- 
ductive gas to the high velocities required. Examples of such studies are 
found in references 1 to 9» Several solutions of the equations for this flow 
are available, such as the one for a constant area and constant temperature 
given in references 1, 2, or 8. Another example is the solution for constant 
pressure given in reference 9 * However, no systematic study of the various 
possible modes of operation of such a channel appears to have been made. The 
results of such a study should be useful for the understanding, evaluation, and 
design of plasma accelerators of this type. 

The purpose of this paper, therefore, is to present a study of a number 
of solutions of the equations of the one -dimensional flow of an inviscid, com- 
pressible, conducting, single-fluid plasma through a channel having mutually 


perpendicular electric-current and constant-magnetic-field vectors, both normal 
to the channel flow direction. For purposes of illustration and comparison, a 
few calculations are presented for most of the modes of operation investigated. 


SYMBOLS 


a = 


All values are given in rationalized mksck system. 

(7 - l)j. 


- 


B0u r 


O' Jr 


A cross-sectional area 

b = aM 0 2 

B magnetic flux density normal to flow, assumed to be constant 

C constant of integration 

c velocity of sound, (yp/p) 1 / 2 

Cp specific heat at constant pressure 

E ratio of effective electric field within plasma to u<^B 

f( ) function 

i,j,k orthogonal unit vectors 

j current density 

M local Mach number of flow, u/c 

p local pressure 

p = P/Po 

t temperature 

T = t/t Q 

u velocity in x-direction 

U = u/u Q 


(U), 


value of U for a stationary value of y 


max 


2 



X 


flow axis in one -dimensional coordinate system (see fig. l) 


y 


nondimen sional length, Bj Q (x - x Q ) /p Q u 0 2 


a 


7-1 

Bcm 0 


7 ratio of specific heats, considered as 1.4 in all computations 

5 increment 

p local density 

P' = P/P 0 

tf electrical conductivity of gas, assumed constant 

Subscript : 

o an initial or reference condition 


THEORY 


A schematic diagram of the system considered is shown in figure 1. It is 
seen that the electric current and the magnetic field are normal to each other 
and to the flow direction along the channel axis. Although these flows may be 
expressed with the Navier-Stokes equations, analytical solutions of those equa- 
tions cannot be obtained; thus, if such solutions are desired, it is necessary 
to make the following simplifications: 


The effects of both viscosity 
and thermal conductivity are consid- 
ered sufficiently small so that they 
may be neglected. The density of the 
gas must be sufficiently large so that 
the mean free path is much smaller than 
any of the dimensions of the channel. 
The number of particles entering and 
leaving the channel is assumed to be 
constant; that is, the percentages of 
ionization and dissociation do not 
change as the temperature and pressure 
change. The actual electric field 
using both normal and axial components 
is assumed to be applied in such a 
manner (ref. 2) that the current den- 
sity is normal to the axis and is 
assumed to be applied through separate. 



Figure 1.- Schematic diagram of channel. 
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individually powered electrodes distributed over the conducting wall. The 
normal component of the applied electric field is the component of the electric 
field used to determine the current density from Ohm's law. The externally 
applied magnetic field is considered to "be constant, normal to both the elec- 
tric field and the flow direction, and much greater than the induced magnetic 
field due to the current through the gas. Because the magnetic Reynolds number 
is assumed to be much less than one, the induced component of the magnetic 
field is neglected. The electrical conductivity is also assumed to be constant 
throughout the channel. Finally, the effect of ion slip is neglected, although 
it is shown in reference 2 that under certain conditions such effects can be 
fairly large. 

The nature of these assumptions, especially those dependent on constant 
temperatures, such as constant ionization, dissociation and electrical conduc- 
tivity, reduces the applicability of the results presented herein. However, if 
greater accuracy of prediction of results for a specific channel is desired, 
the channel could be broken into a number of segments over which the variation 
of all the parameters involved is sufficiently small so that the results of 
this analysis are applicable to the segment under consideration. The condi- 
tions at the end of the segment may then be used to compute the new initial 
conditions for the next segment. By continuation of this process, a solution 
could be built up for a specific channel which would show the effects of 
changing the magnetic field, the conductivity, the percentage of ionization and 
dissociation, and the ratio of specific heats. Another limitation is seen in 
the fact that if the area changes are too great, serious departures from the 
one -dimensional treatment given herein will occur. For example, in the case 
where the flow is supersonic, the resultant wave structure, especially if the 
waves are shocks, can cause such serious departures of the properties of the 
actual flow from the theoretical flow that the results of an analysis would be 
made inapplicable. Thus, care must be taken in applying the results of the 
variable-area studies to any channel design to use values of the parameters B, 
M q , and p 0 u Q 2 for which the rate of change of area with respect to the actual 
channel length is sufficiently small to reduce the two- or three-dimensional 
effects to values that are small compared with the one -dimensional effects. 

With proper consideration of these assumptions, the equations for the 
quasi one -dimensional flow of an inviscid gas through crossed electric and 
magnetic fields are 


^ + = 0 
p u A 

which expresses the continuity of the flow, 

pu du + dp = jB dx 

which expresses the change of momentum. 


pu d^Cpt^ - u dp = 


j 2 dx 


(la) 


(lb) 


(lc) 


4 


which represents the heat added to the stream, and 


. S£ + “ (id) 

p p t 

which is the equation of state. The additional relations derived from the 
definition of M and relations (la) and (id) 

2 dM _ du dp dA _ 2 du dt /, \ 

M u“p“A u " t v ' 


are convenient for determining the Mach number M. It is convenient to elimi- 
nate du and dt from equations (lb) and (1c) with the use of the values of 
du and dt from equations (la) and (id). If these eliminations are made and 
c2 (the square of the velocity of sound) is substituted for 7p/p, it follows 
that 

-P U - ^£ + dp = JB dx + pu 2 (2a) 

and 

dp + u dp - - 11 J 2 ax (2b) 

P $ 


These equations may be solved for dp and dp to obtain differential equa- 
tions in dp and dp. Corresponding equations for du, dt, and dM are 
obtained by use of equations (la), (id), and (le). 


If proper manipulations are made on these differential equations and the 
quantities U = P = P* = T = ^, and a = 5“^ are introduced, 

the following equations result: 


dU ... 
U dx 



M%JA 

Po^o^Aq 


<U - aj) + 


dP 

P dx 


7M 2 

r b.ia 

M 2 - 1 

Po u o 2ij2a o 


(U - aM 2 j) 



(5a) 


(5b) 


dT = (7 - l)M 2 f BjA / 

T dx M 2 - 1 [p 0 u 0 2 U 2 A 0 \ 



( 5 c) 
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dp' 


(5d) 


p' dx 


m 2 r 

= -iL; 


BjA 


PoV^ 


-(U - aj) + 


(3A 

A dx 


3 


dM 

M dx 


f M%jA 1^ ( 7 M 2 + l)aj 


2 + (7 - l)M 2 

Wuo^Aol 7 + 1 J 

T 

7+1 


dA 


A dxf 


(3e) 


The differential relations (la), (id), and (le) may be integrated to give 

p'U— = 1 (4a) 

A o 


P = p'T 

M 2 _ u£ _ U 2 P' _ UA o 

M 2 T P PA 
o 

With the use of these relations only equations (3a) and (5b) need to be inte- 
grated once the area distribution A is known. Inspection of equations (3a) 
to 3(e) shows that there are two more dependent variables (that is, U, P, 
p', T, M, A, and j) than there are equations; hence, it is necessary to 
specify any two of these variables before solutions can be obtained. For this 
analysis either the variation or the constancy of two of them will be specified. 

Several important features should be noted concerning the limits of integra- 
tion and methods of computing any of the solutions considered in this analysis. 
Since the flow must always start from the initial conditions and proceed to the 
local condition under consideration, the limits of integration must also extend 
from the initial value of 1 to the local value; or if the equation contains a 
constant of integration, that equation must be evaluated at the initial flow 
values. Also since the flow cannot reverse upon itself, all computations must 
start from small increments of the initial values and then follow the direction 
required to render x positive. This procedure is necessary since in many of 
the equations it will be convenient to treat U or P as the independent 
variable and to solve for x. Thus, if U is considered to be 1 + S where 
8 is small compared with 1 , the sign of 8 must be chosen to render x 
positive, and then the computations may continue with U decreasing or 
increasing as indicated by the sign of the term 8 . 

A general note should also be made concerning the singularities which 
appear in many of the integrals required for the evaluation of x. In all such 
integrals the value of the integral approaches infinity as the variable of 
integration approaches the singularity. Since x cannot proceed beyond an 
infinite value, the integration need not be carried beyond the first singular- 
ity. Care must be taken therefore that the integrand is not infinite for any 


(»*) 

(4c) 
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value of the variable of integration falling between the considered local value 
and the initial value. 


VARIABLE -AREA CHARNELS 


Constant Velocity 


The velocity in a channel may be held constant by allowing the area 
to vary such that dU is always equal to zero. Such an area variation 
is attained by setting the term within the brackets of equation ( 3 a.) equal 
to zero. 


dA _ M g ABj(l - aj) ^ 
A Po u o 2a o 


( 5 ) 


wherein the constant value U = 1 is substituted for U. If equation (5) and 

2 

the relation M g = ^ from equation (^-c) are substituted into equa- 

PA 

tion ( 3 b), it follows that 


dP 

P 


TM^ 

P 


BJ 

Po u o 2 


dx 


(6) 


Since the flow in any channel must necessarily start at the reference condi- 
tions, the limits of integration must cover the region from P = 1 to P and 
from x = x Q to x. With the velocity U constant, j cannot be a function 
of U; however, the pressure integral can be readily evaluated if j is chosen 
as a function of P such as j = j Q f(P) where f(l) = 1 . The solution for 
the condition that U = 1 becomes 


1 r P dP 

7Mq2 J 1 f (P) 


jp 

p u 2 

M o o 


(* 


x o) = y 


( 7 ) 


where y is a nondimensional quantity proportional to the channel length. The 
parameter y will be used throughout this analysis to compare channel flows 
under various conditions. 

If equation (5) is expressed in terms of dP/P rather than in terms of 
dx, the result is 

dA _ J- ~ a j dp / n\ 

A " 7 P K } 

P, equation (8) may be integrated and the 


Again, if j is a known function of 
limits evaluated to give 
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(9) 



f*\} - ¥ 


Thus A can be expressed as a function of P and therefore of x. 


Constant Pressure 


General solution . - For the general solution the pressure will be assumed 
to be constant (P = l) throughout the channel. This condition appears to be 
important since this is the area distribution that the flow will take if the 
channel walls are free fluid surfaces. Thus, these walls are naturally- 
occurring walls and no artificial restraint need be used to cause the flow to 
conform to a given area distribution. In all the other cases, the flow must be 
constrained to a given tailored wall distribution. The pressure ratio P may 


be rendered constant by setting 
stitution of = 


_ M q "UA 0 



in equation (3b). 


Then with the sub- 


BjA (a ~ ) 

P 0 Uo 2 A 0 U 


dx 


( 10 ) 


With the substitution of equation (10) into equation (3a) it follows that 


dU 


BJA 


p o u o “o 


2a, 


dx 


(H) 


The elimination of 


-®iA_dx 

Po u o 2 Aq 


from equations (10) and (ll) results in 


dA _ dU 

A - aMo2 Ao j U 

This equation is a first-order differential equation of the form 

d4 + A = ^AqJ 
dU U U 


( 12 ) 


(13) 


Its integral is given in most standard differential equation texts (ref. 10, 
for example) and is 


_A_ 

Aq 



( 1 *) 
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The velocity is now given by equation (ll) 


as 


au 


B 

Po u o 2 


^ M o 2 

U 



( 15 ) 


In the three subsequent sections j will be taken as the following functions 
of U: j = j Q , j = j D U, and j = Bcru 0 (E - U) where E is constant. The 

current distribution j = Bcru 0 (E - U) is the simplified form of Ohm's law 

expressed in terms of the nondimensional parameters of this paper. If E, 
which is proportional to the electric field within the plasma, is held con- 
stant along the channel, it may then be considered that the electric field 
within the plasma is held constant. Also the problems are treated wherein the 
current density varies inversely as the cross-sectional area of the channel. 

Constant current density, j = j . - The channel which has a constant cur- 
rent density, j = j 0 , is one of the more practical of the channels investi- 
gated in this paper because of the ease of approximating constant current den- 
sity. One relatively simple method of approximation of constant current 
density would be the location of a large number of electrodes, or arcs, over 
the channel surface such that the area allotted to each electrode is constant. 
Then, of course, it is necessary only to adjust the driving potential on each 
electrode so that the currents are the same. If j = j D and if the term 

0 tM o 2 j o is set equal to b, the evaluation of the integrals in equation (l4) 
results in the following area distribution: 

A_ = b(U - 1) + 1 06) 

U 


The substitution of j = j Q and b = aM 0 2 j 0 into equation (15) and the subse- 
quent integration result in the following nondimensional velocity distribution: 


y = 




X* 


PnU, 


o“o 




— ^log[b(U - 1) + 


(17) 


Since P = 1, the relations between the remaining dependent variables from 
relations (4c), (4b), and (4a) follow: 


M2 _ UA 0 
M 0 2 A 





J 


(18) 
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Several velocity, area ratio, and Mach number distributions computed from 
equations ( 17 ), (l6), and (l8), respectively, are presented in figures 2, 3> and 

4. It is interesting to note that 
the velocity increases regardless of 
the value of b, which includes ini- 
tial Mach number, magnetic field, 
conductivity, and initial velocity. 
(See fig. 2.) This mode of operation 
is therefore indicated if velocity 
acceleration is desired in the chan- 
nel. It is also observed that if b 
is made sma.ll, that is, less than 
about 0, then the accelerator becomes 
insensitive to changes in b. How- 
ever, if area changes are considered 
to be undesirable, operation in the 
region of b = 1 for which the area 
is constant would seem more desirable, 
especially if a smooth and -uniform 
flow is desired. The areas are also 
seen to approach the constant b as 
y or U increases. (See fig. 3* ) 

In equation (l6) if U becomes large, 
A/A 0 = b. 

Examination of figure 4- shows 
that the Mach number increases. How- 
ever, the trend for large values of b 
(see b = 10 curve in fig. 4) indi- 
cates a reduction in Mach number which 
will be recovered later. This opera- 
tion should also be desirable for a high Mach number output since Mach number 
ratio gains of 2 to 3 are reasonable. For example, if the input Mach number is 4 
or 5 , output Mach numbers of an order of 10 to 20 should be reasonably attain- 
able, especially if the conductivity and initial velocity are sufficiently large, 
so that b is of the order of 0.1 or thereabouts. Also the Mach number curves 
indicate that Mach number can be controlled by changing the value of b; however, 
such changes would also change the final velocity, and only a detailed study 
would show whether such an operation is feasible. The reduction in Mach number 
for the b = 10 curve occurs because the increase in temperature which is pro- 
portional to UA/Ao is sufficiently greater than the increases in U and A 
to lower the value of the Mach number. 

Current density proportional to U, j = j Q U. - The mode of operation when 

the current density is proportional to U, j = j Q U, is included because of its 
mathematical interest rather than because it represents a practical channel. 
However, as is shown under the special cases, this mode is required in a 
constant-area, constant Mach number channel. This mode is impractical because 
it requires velocity measurements and tailoring of the electrode current readings 
to those measurements and because, as is shown by the analysis, the velocity and 
current approach infinity within a finite accelerator distance, a condition 
impossible to meet. 



Figure 2.- Velocity ratio distribution in a 
magnetohydrodynamic channel with P, B_, 
a, and j held constant. 
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The flow conditions are determined by the substitution of j = j Q U into 
equations (l^) and (15) ■ Evaluation of appropriate limits and constants of 
integration results in the following distribution for a/A 0 and U: 



If b < 2, 



x 0 ) 


2 

\|b /2 - b 





(20a) 


if b = 2, 



(20b) 


and if b > 2, 



Equation (l8) again follows for the other dependent variables. 


(20c) 



0 .4 .0 1.2 1.6 2.0 2.4 

y 


Figure 5-- Velocity ratio distribution 
in a magnetohydrodynamic channel with 
P, B, a, and j/U held constant. 


Computed values of U, A/A 0 , 
and M/Mo f° r various values of b 
are presented in figures 5* 6, and 7> 
respectively. The surprising point 
seen in the velocity curves (fig. 5) 
is the fact that for values of b > 0, 
the velocities approach infinity at 
some definite value of y. This phe- 
nomenon appears to be reasonable when 
it is recalled that the current den- 
sity increases with U; hence, the 
acceleration forces also increase. 
Thus, the two combine to produce an 
infinite velocity. This phenomenon 
is also easily seen in equations (20a) 
to (20c) by letting U become very 
large; then y approaches 



various ranges of b. 
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Figure 6.- Area ratio distribution in a magnetohy- 
drodynamic channel with P, B, a, and j/U held 
constant . 



Figure 7- - Mach number ratio distribution in a 
magnetohydrodynamic channel with P, B, a , 
and j/u held constant. 
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Correspondingly, it is seen that the area ratio (fig. 6) also approaches 
infinity. Since the area and velocity both become infinite for this current- 
density distribution at some given length, any channel involving this mode of 
operation must necessarily be ended before that length is attained. 

This mode of operation will increase the Mach number ratio (fig. 7) provided 
b is kept less than 2. Greater values of b will reduce the Mach number ratio. 

Constant electric field, j = Bctuo(E - U) . - Operation of a channel with a 
constant potential is often proposed because in the ideal channel only the 
electrode voltages need to be held constant. This mode of operation has the 
advantage at high Mach numbers that as M increases, the heating rate j 2 /cr 
decreases. This property is desirable because at high Mach numbers a high 
heating rate will rapidly increase the entropy. In the application of these 
constant electric field equations, the actual electric field within the plasma 
should be used rather than the electric field given by the potential between 
the electrodes, because of sheaths. 

The equations representing the effects of a constant electric field are 
obtained with the substitution of j = Bauo(E - U) into equations (l4) and 
(15) • These substitutions result in the following relations: 

^ = ~ - D(2E - U - 1) +0 (21) 


y = 


B Jp 

Po u o 2 


(* 


x o) 


. E ~ 1 2 r o f- 2E lo S 

(7 - 1)Mo 2 ^L 


E - U 
E - 1 


+ (E + Q) log 


E + Q - U 
E + Q - 1 


where 


+ (E - Q)log 


E - Q - U 
E - Q - 1 


Q 2 = (E - l) 2 + 


2 

(7 - D^ 2 


( 22 ) 


(23) 


and 


= BoUq(E - 1) 


(24) 


It is observed that these equations are functions of two parameters E and 
M 0 rather than of the single parameter b as in the previous cases. Thus, 
there will be a family of constant E curves for each initial Mq or a family 
of constant Mo curves for each initial E. Plots, of a fe\f computed values of 
U, k/ko, and Mach number obtained from equations (22), (21), and (l8) are 
shown in figures 8, 9, and 10, respectively. The velocity ratios (fig. 8) are 
seen to approach E, as may be expected, since the current density is zero if 
U = E and no further acceleration can occur. Thus, the constant-pressure 
mode with a constant E is l imi ted in its acceleration. This is not, however, 
serious since E can readily be made large. 
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Figure 8.- Velocity ratio distribution in a 

magnetohydrodynamic channel with Pj B, a, 
and E held constant. 



Figure 9 . - Area ratio distribution in a 

magnetohydrodynamic channel with P, B, 
a, and E held constant. 


The final velocity ratio can "be 
shown from equation (22) to "be inde- 
pendent of the initial Mach number. 

The initial Mach number does, however, 
affect the distance in which the veloc- 
ity ratio approaches its limiting value. 
These curves indicate that a high ini- 
tial Mach n umb er should be used if rapid 
acceleration of the gas is required. 

It is interesting to note that the 
conductivity appears only in the term 
j = Bau 0 (E - l) of the scale factor 

Bj 0 /p 0 u 0 2 which converts y into an 
actual length. 

Readily traceable trends are not 
so easily seen in the area distribution 
and the Mach number curves (figs. 9 and 
10, respectively) except that both of 
them tend to approach constant values 
as y increases. The Mach number is 
seen to increase for all values of E. 



Figure 10.- Mach number distribution in a mag- 
netohydrodynamic channel with P, B, a, 
and E held constant. 
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However, for the condition that Mq = 2 and E = 10, the Mach number first 
decreases (a condition indicating a rapid heating of the plasma) before it 
increases to a value greater than the initial Mach number. 

Current density inversely proportional to the area, jA = JqAq.- The opera- 
tion of a channel with the current density inversely proportional to the cross- 
sectional area has an interesting feature. If the channel height is constant, 
the current through any section of the channel of length Sx is equal to the 
current through any other section of the same length. Such a channel could be 
called a constant -current channel. This current distribution is approximately 
attained by equally spacing the electrodes or sets of electrodes along the 
channel axis and making the currents in each set equal to the current in every 
other set. 

Since j is a function of A rather than of U, equations (l4) and ( 15 ) 

are not directly applicable. However, the substitution of j = 
tion (ll) results in 




into equa- 


dU 


PoV 


dx 


ay 


(25) 


which may be immediately integrated to obtain, after evaluation of the limits. 



Figure 11.- Velocity ratio distribution of a 
magnetohydrodynamic channel with P, B, 
a, and JA held constant. 


y = U - l (2 6 ) 

or 

U = 1 + y 


Also the substitution of JA = j 0 A 0 
into equation (12) results in 


A 

A 2 - bA 0 2 



(27) 


which upon integration and evalua- 
tion of limits and after substitu- 
tion of equation ( 26 ) becomes 

= [b + (l-b)(l+y)- 2 ] l/2 ( 28 ) 
A o 


Typical computed values of 
the velocity ratios, area ratios, 
and Mach number ratios are pre- 
sented in figures 11, 12, and 13 , 
respectively. As seen from equa- 
tion (26) or figure 11, U is a 
linear function of y with no 
other parameters involved, except. 
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Po u o 

of course, the scaling factor -— 5 — 

■“Jo 

which converts y into actual length. 

The area ratio (fig. 12) does involve 
the parameter b, the ratio increasing or 
decreasing depending on whether b is 
greater than or less than 1. Also the 
area ratio appears to approach a constant 
value as y becomes large. This value 
is seen to be equal to b -^-/ 2 by letting 
y be large in equation (28). This mode 
of operation increases the Mach number 
for all the values of b considered, 
that is, 10 or less. It can be seen 
however that regardless of the value of 
b, the Mach number ratio will eventually 
increase provided y is sufficiently 
large. This is true because by equa- 
tion (l 8 ) M = Mq^UAq/A; therefore, A 



Figure 12. - Area ratio distribution in a 
magnetohydrodynamic channel with P, 
B, cr, and jA held constant. 


eventually becomes constant but U continues increasing and M must eventually 
vary as Ul/ 2 . Thus, M/M q must finally become greater than 1. An accelerator 
operating under this mode should increase the Mach number substantially if b 
is rendered small by a large value of the conductivity a. 



Figure 13- - Mach number ratio distribution in a 
magnetohydrodynamic channel with P, B, a, 
and JA held constant. 
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Constant Density Solutions 


General solution .- The density ratio p' rather than the pressure ratio 
P may he held constant. This mode of solution may he obtained hy determining 
the value of dA/dx from equation (3d) with the condition that dp' equals 
zero. However, the area ratio for constant density is more readily determined 
from equation (4a) . Since the constant density i's expressed hy letting p = 1, 
equation (la) becomes 


dA = dU 
A “ ~ U 


(29) 


which upon integration becomes 


A_ _ 1 
A 0 U 


(30) 


The use of equation (29) in equations (3a) and (3b) will give the equations in 
dU and dP which express the condition that p is held constant. After 
these substitutions and appropriate reductions are made, it is found that 


II 2 

U - aj 


dU 


Bj 

Po u o 2 


dx 


and 


(3D 


U 


o,7M 0 2 j 


dP 


= 

Po u o £ 


dx 


(32) 


The elimination of 


Bj 

Po u o 2 


dx 


from equations ( 31 ) and ( 32 ) results in 


7M 0 2 


dP = 


oUj 

U - aj 


dU 


(33) 


It is seen by inspection that the variables U and P are separated if j is 
a function of TJ and hence P can be expressed as an integral. 

Constant current density, j = j Q . - In this section, the current density 

is assumed to be constant along the length of the channel. The substitution of 
j = j 0 into equations ( 31 ) and (33) results in 
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dx = dy 


and 


U 2 

U - a 


dU 


Bjp 

Po U o 2 


(54) 


1 

7M 0 2 


dP = 


ajpU 
U - aj Q 


dU 


aU 

U - a 


dU 


(55) 


where a = <xj 0 = 


(y - i)3 c 

Bcru„ 


If equations (54) and (55) are integrated and the 


constants of integration evaluated from the appropriate reference values, then 


y = u2 ~ ^ + a(U - 1) + a 2 log ^ a 

cL - a 

and 


(56) 


— - = a(U - l) + a^ log — ~ a 

7M q 2 1 " a 


(57) 


(7 - l)j 

where a = aj 0 = — . It is seen that U is independent of M Q , whereas 

BctUq 

P depends on Mq as well as on a and U. 


The other variables, M, T, and A, are obtained from equations (4c), 
(4b), and ( 30 ), respectively, as 


2 _ Mp IF 

P 

( 58 a) 

T = P 

( 58 b) 

> 

11 

( 58 c) 


These equations show that the Mach number and the temperature ratio of the 
flow in a magnetohycLrodynamic channel in which the density and the current den- 
sity are held constant are dependent on the initial Mach number, whereas the 
area ratio is not. 
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P-1 


as functions of a 


Several curves showing values of U, A/A 0 , and. 

7Mq 2 

and y are presented in figures l4, 15, and 1 6 , respectively. Examination of 
the curves (fig. l4) presenting the velocity ratio U shows two important fea- 
tures. First, that acceleration occurs only if a is less than 1 in contrast 
to the constant-pressure mode of operation in which acceleration occurs for all 
values of b or aMo 2 . Second, it is seen that if a is greater than 1, the 
flow velocity is reduced to zero. Such reduction involves, of course, the 
filling of an infinite sink since the area becomes infinite when the velocity 
is reduced to zero. Thus the flow is literally stopped in this mode of opera- 
tion if a is greater than 1. Also, it is observed that the velocity ratio 
becomes a very sensitive function of y whenever a has values close to 1 as 
indicated by the large changes of U with respect to y when a changes from 
O .99 to 1.01. 

The area ratio curves (fig. 15) are obtained from the reciprocal of the 
velocity ratios (fig. 1^) and show that for a > 1, the area ratio approaches 
infinity as U approaches zero and y remains finite. In figure 16, which 

P 1 

presents the pressure ratio function , it is seen that P and hence T 

7Mo 2 

remain finite as U approaches zero when a is greater than 1. Hence, if a 
is greater than 1, this mode of operation not only reduces the velocity to zero 
but also restricts the pressure and temperature ratios that can be obtained. 
However, the values of U, P, and T can all be increased if a and B are 
sufficiently large so that a is less than 1. 



Figure Pt-.- Velocity ratio distribution in a mag- 
netohydrodynamlc channel with p ' , B, a, and 
j held constant. 
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A. 

Ao 


Figure 16.- Distribution of pressure 
function in a magnetohydrodynamic 
channel with p' B, a, and j 
held constant. 


Figure 15. - Area ratio distribution 
in a magnetohydrodynamic channel 
with p', B, a, and j held 
constant. 




Current density proportional to velocity, j = j Q U.- In this section the 
current density will he considered to vary directly with U (that is, j = j 0 u). 
Moreover, because in a constant-density channel, U varies as A 0 /A, the rela- 
tion jA = JoA 0 is also satisfied. 

The functions U and P are obtained with the substitution of j = j D U 
into equations (31) and (33)- This substitution results in 

■ dU = ay (39) 

1 - a 

and 

= - a — p dU (40) 

7M q 2 1 - a 
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in 


The integration of equation (39) and. the evaluation of the constant result 



or 


U = 1 + (l - a)y 


(41) 


The integration and evaluation of equation (4o) with the substitution of equa- 
tion (4l) result in 


P - 1 _ ay {g + (1 - a)yl 

/Mq 2 2 


(42) 


Equations (^l) and (42) along with equations ( 38 a) to ( 38 c) fully express the 
flow properties of this particular mode of operation of a magnet ohydrodynamic 
channel, that is, with p', B, and cr constant. 

Several curves are shown in figure 17 which presents U as a function of 
y for various values of a. Again, U "becomes zero within a finite distance 
for values of a greater than 1 "but accelerates without limit for values of a 
less than 1. Correspondingly, the values of the area ratio A/A 0 which are 
not shown either increase to infinity or asymptotically approach zero as y 
"becomes infinite, depending upon whether a is greater than or less than zero. 



(a) a ^ 1. (t>) a i 1.1. 

Figure 17 .- Velocity ratio distribution in a magnetohydrodynamic accelerator 
with p 1 , B, a, and j/u held constant. 
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Values of the pressure function for a less than 1 are shown as functions 
of y in figure 18(a) and for a greater than 1 in figure l8(h) . If a is 
less than 1, it is seen that the pressures increase indefinitely. If a is 

greater than 1, then, as in the previous case, the pressure function ^ 

7Mq 2 

and, correspondingly, P and T rise to a limited but finite value. There- 
fore, this mode of operation has several characteristics similar to the con- 
stant current density mode in that the value of a strongly determines the 
nature of the operation of the channel. 



Figure 18. - Pressure ratio function in a mag- 
netohydrodynamic channel with p', B, a, 
and j/U held constant. 
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(b) a S 1. 

Figure 18. - Concluded. 

Current density inverse ly proportional to U, jU = j Q . - In this section, 
the current density is considered to vary inversely with respect to U. 

Such a variation may he expressed as j = —2.. The values of U and P for 

U Jo 

this current distribution are determined with the substitution of j = — 

into equations (31) and (33). After making the indicated substitution, equa- 
tions ( 31 ) and (33) become, respectively. 


— dU = dy (43) 

Ij2 - a 

and 

— dP = — — — dU (MO 

7M q 2 U 2 - a 


After equations (4-3) and (44) are integrated and the constants evaluated, 
they become, respectively. 


y = + a(U - 1) 

3 


+ 



2 


log 


(u - ^)(l +y/i) 
(u + ^)(l - y/l) 


(45) 
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and 


P^-l.= |lo g n!_^_a (46) 

7 M Q 2 2 1 - a 


Equations ( 58 a), ( 58 b) , and ( 38 c) may again be used to determine M, T, and 
A/Aq. 


Values of U and — as functions of y are presented in figures 19 

7M 0 2 

and 20. Examination of figure 19 shows that for this mode of operation if a 



is less than 1 , then 


y = H!_i_ 

* 3 


and 


the curve for a = 0 is an excellent 
approximation of the velocity distri- 
bution. However, when a becomes 
larger than 1 , an almost discon- 
tinuous change occurs in which y 


changes from the 


u3 - 1 
3 


distribu- 


tion to a distribution in which U 
quickly becomes zero, for example, 
see the a = 1.0001 curve; or if 
a = 1.225, then U = 0 at y = 0.1^. 
The areas become infinite just as in 
the previous cases. Perhaps, the 
interesting point is the fashion in 
which the velocity curves approach 1 
as a approaches 1. This phenomenon 
is indicated at the base of the 
a = 0.9998 curve. This curve indi- 
cates that as a becomes still closer 
to 1, the U curve will fall on 
the value U = 1 for larger and 
larger values of y until finally it 
breaks away and becomes approximately 
parallel to the a = 0 curve. 


Figure 19. - Velocity ratio distribution 
in a magnetohydrodynamic channel with 
p', B, a, and jU held constant. 

increase for all values of a less than 


The pressure function curves for 
a < 1 are shown in figure 20. These 
curves show that the pressures 
1 and greater than 0. At a = 1, the 


pressure function is expressed by 
stitution of j = — and U = 1 

u TT 


r ~ x = y. This may he shown with the sub- 
7Mq 2 

into equation (32). If a is greater than 1 , 


the pressure curves are terminated at a finite value of y. 
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Figure 20.- Pressure distribution function in a mag- 
netohydrodynamie channel with p', B, a, and jU 
held constant. 


Constant electric field, j = Bcru 0 (E - U) . - The current density which 
occurs if the electric field within the plasma is constant has "been shown to 
be j = Bctu 0 (E - U) . The use of this function in equations (jl) and (33) 
results in 


U 2 dU 


= B 2 o 


dx = 


BJ, 


O - (7 - 1)E](U - E) Po^o p 0 u q 2 (E " 1) 


dx = - 


By 

E - 1 


and 


(vr) 


dP 

7Mo 2 


(7 - 1 )(E - U)U 
7 U - (7 - l)E 


dU 


m 


After integration and evaluation of the constants, equations (47) and (48) 
become, respectively. 
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y = (E - 1) 


'il - 1 ) g E log w --Jl - - E log - 2^ 1 

7 2 y - (7 - 1)E 1 - E 7 


3 


and 


P - 

7Mq : 


-i = I ~ 1 T 

- 27 


1 - U' 


.2 + 2E(U - 1) + 2(7 - 1)E 2 7 U - (/ - 1)E 


7 - (7 - 1 )E 


E 


(^9) 


(50) 


It is interesting to note from equation (49) that U is constant for two 
values of E; that is, E = 1 and E = — — _ (e = 5.5 for 7 = 1.4). It can 
also he shown that U approaches E asymptotically as y becomes infinite, 
provided 1 < E < — 2L_. These observations are verified upon examination of 


figure 21, which presents values 
of U plotted as a function of 
y for different values of the 
parameter E. If E is less 


than 



but greater than 1 


= 3-5 if 7 = 1.4^, the 

U-distribution approaches E as 
y becomes large. However, it 
is also observed that as E 
closely approaches 5*5* the veloc- 
ity U follows the U = 1 curve 
for a considerable length before 
it leaves the curve and approaches 
the E = 5.5 value. If E = 3-5, 
the velocity becomes constant at 
U = 1 as may be seen by examina- 
tion of equation (47). In this 


equation, if E = — - , dU must 

7-1 

be equal to zero. Thus, since dU 
equals zero, U must equal 1. 

Eor larger values of E, dXJ must 
be negative; thus, as y increases, 
U must decrease as is indicated by 
the curves. Again, infinite cross- 
sectional area is required to meet 
continuity conditions. The inter- 
esting point here is that if E 
is made large, the flow is stopped 
rather than accelerated to the 

value of E or even — as 

7 

might be expected from equation (49) . 



(a) e < 3 . 5 . 



(b) e > 3.5. 

Figure 21.- Velocity ratio distribution in a 
magnetohydrodynamic channel with p ' , B, 
and a held constant and j = B<jU 0 (E - U) . 
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The corresponding pressure function 


P - 1 


is presented in figure 22. 


7M 0 2 

These curves show that the pressure function approaches a limiting value for 
values of E between 1 and 3*5* If E is greater than 3-5* "the pressure 
function terminates, as might be expected, at the same value of y at which 

P - 1 

U becomes equal to zero. At E = 3-5> the value of = is determined by 

7M 0 2 

substitution of U = 1, j = j 0 = Bcruo(E - l), and E = — I — or 3-5 into 

7-1 

equation (32). The result after integration and evaluation of constants is 



(51) 



Figure 22. - Pressure function distribution in a 
magnetohydrodynainic channel with p', B, and 
a held constant and j = BffU 0 (E - U) . 


Constant Temperature 

•The operation of a channel at constant temperature is often proposed (see 
refs. 1 and 2 ) because, among other reasons, the gas properties such as con- 
ductivity, ionization, dissociation, and the specific-heat ratio 7 remain 
constant. This mode of operation is usually obtained by properly varying the 
current density in a constant-area channel. It is, however, possible to vary 
the area in such a fashion to insure constant temperature. 
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This constant -temperature area variation is obtained by setting 


1 dA - BjA L 
A ^ Po^o^Ao V 



( 52 ) 


which, upon substitution into equation (3c) renders dT equal to zero; and 
hence,, T becomes a constant. Since T is defined as a temperature ratio, 
its value must equal 1. Also since T = 1, M = M 0 U, from equation (4c) . With 
the substitution of M = M 0 U, equation (52) becomes 

1 dA _ BjA 
A 4x1 PoUo^Ao 

If equation (53) is used in equation (3a), that equation becomes 


7 M Q 2 U 2 aj - (7 - 1)U - 




(53) 


The elimination 
lowing relation 


dU 

U dx 


of 


BjA 


Poiio^A 

between U 


BjA 

Pq u o^^A 0 


(7 - 1)U + aj 

7-1 


(54) 


— from equations (53) and (54) results in the fol- 
o 

and A .’ 


dA = (7 M q 2 U 2 - l)aj - (7 - 1)U dU 
A (7 - 1)U + qj U 


(55) 


If j is a function of U, the variables in this equation are separated so 
that the solution of equation (55) is reduced to a quadrature. Once A/A 0 is 
obtained, its value may be substituted into equation (54) which may then be 
integrated to obtain both A/A 0 and U as functions of x. Unfortunately, 
the simple current-velocity relations used previously do not give functions of 
A/A 0 which, when substituted into equation (54), will render that equation 
integrable with elementary functions . 


Constant Mach Number 


General solution .- The general solution for M = Mq is obtained in a 
similar fashion. If cLA/A is given by 


1 dA _ 
A dx 


M 0 2BjA 

(7 + 1)U - ( 7 M 0 2 + l)qj 

PoUq^Ao 

L 2 + (7 - 1 )Mq 2 


(56) 
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the substitution of this equation into equation (3e) will make dM = 0, and 
hence M is a constant. The substitution of equation (5 6 ) and M = M 0 into 
equation (3a) gives the following equation for dU: 


dU 

U dx 


Mo 2 BJA 

(7 - 1 )U + aj 

Po^o^Aq 

2 + (7 - a)M 0 2 


(57) 


The elimination of 
results in 


Mp^BjA 

P 0 u 0 2 u2Ao[2 + (7 - 1)M q 2] 


from equations ( 5 6 ) and (57) 


dA _ (y + i)u - (yM Q 2 + l)a-j au 
A " (7 - 1)U + aj U 


(58) 


Again, if j is a function of U, equation ( 58 ) is separated and the solution 
is reduced to a quadrature. Again U is determined by integrating equa- 
tion (57) after substituting the value of A/A 0 as determined by the integra- 
tion of equation ( 58 ) as follows: 


U dU 


M d 2 B dx 


J [(7 ~ 1)U + aj] e 


T (7+ l)U-(7M 0 2 +l)qj 
(7-l)U+aj 


dU 

U 


+C 


p 0 u 0 2 [2 + (7 - 1)M 0 2] 


(59) 


Current density proportional to U, j = j Q U. - Of the five functions of 
U that have been used to express j, only one, current density proportional 
to U, that is, 

J = Jo 11 (80) 


will render equation (59) readily integrable. The substitution of this func- 
tion into equation ( 58 ) gives 


dA _ y + 1 ~ (y M o 2 + 1 ) a dU 
A ~ 7 - 1 + a U 

A similar substitution into equation (57) gives 


dU 

U 


= M q 2 


7 - 1 + a _A_ 

_2 + (7 - l)M 0 2 A q 



(61) 


(62) 
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If equation ( 6 l) is integrated and the constant evaluated, it becomes 

7 +l-(yM Q 2 +l)a 


A _ 


= U 


7 - 1 +a 


(65) 


If equation (65) is substituted into equation (62), that equation becomes 


7+l-(7Mo 2 +l)a n 


u 


7 - 1 +a 


M 0 2 (7 - 1 + a) 
dU = - 2-11 . — dy 


2 + (7 - 1 )M 0 2 

After equation (64) is integrated and its constant evaluated, it becomes 


(64) 



2 + (7 ■ 

- 1 )M 0 2 L 

Mo 2 ! 

[V + 1-1 

JyM 0 2 + l) 

,a]\ 


7 +l-( 7 M Q 2 +l)a 

7 -l+a 


- 1 = y 


(65) 


or 


7 - 1 +a 


U 


= [± + M o g [r + 1 - (^Mp 2 + l)a] y . 


7+1- ^7M 0 2 +l) a 

2 + (7 - 1)M 0 2 

Substitution of equation (66) into equation (65) results in 


(66) 


JL = 


2 + (7 - 1)M 2 


A ° 2 + (7 - 1)Mq 2 + M 0 2 [ 7 + 1 - ( 7 Mo 2 + l)9y 


(67) 


The variables P, T, and p' may be determined from equations (4a), (4b), and 
(4c) by letting M = Mq. They are 


P = U— 


A 

A 


T = U c 


p. = 

AU 


(68) 

(69) 

(70) 


Examination of equations ( 66 ) and (67) shows that U and A/A 0 are 
dependent on both a and Mq 2 . It is also seen that equation ( 66 ) becomes 


31 



indeterminate at a = . This indetermination is resolved with the use 

7Mq 2 + 1 

of the given value of a and equations (62) and ( 63 ) • The solution of these, 
equations then gives 


7Mq 2 

U = e 7M ° 2+1 


y 


(71) 


and A becomes constant. 


CONSTANT -AREA CHANNELS 


General Solution 

The equations which express the properties of the flow in a constant-area 

magnetohydrodynamic channel are obtained by making equal to zero in 

A dx 

equations (3a) to (3e); this procedure renders A constant and makes A/A q 
equal to 1. With these substitutions, the flow equations for a constant-area 
channel become 


dU = dp' 

U dx p' dx 


M 2 Bj 
M 2 - 1 p 0 u 0 2 U 2 


(U - aj) 


( 72 a) 


dP 

P dx 



Bj 

Po u o 2u2 


(U - oM 2 j) 


(72b) 


dT 

T dx 


dM 
M dx 


M 2 

M 2 - 1 


Bj 

Po u o 2u2 



M 2 Bj ["/ 

2 (M 2 - l) p 0 u 0 2u 2 L 


- 1)U - ( 7 M 2 - l)ajj 
7 + l)U - ( 7 M 2 + l)aj 


and equations (4a), (4b), and (4c) become 


P' 


1 

U 


T = UP 



( 72 c) 


( 72 d) 


(73a) 

(73b) 

(73c) 
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If j is a function of either U or P, the integrals of equations (72a) 
and (72b) may he written as quadratures. These quadratures are obtained by 

first eliminating M 2 from both equations with the substitution of M 2 = . 

This substitution results in the following relations for dU and dP; 


UMq 2 - P M 2 ^ 

- 7 T— r ' du = j2 -4 to 

U - aj PoUq 2 


™q 2 - P ^ = 7M 0 2Bj 
P - aMo 2 J 


dx 


Po u o £ 


r\ 

If the term — 4 dx is eliminated from equations (7^0 and (75) > then 

Po u o 2 


h. - °Q) 


dP = y F ~ ^o 
dU U - aj 

This equation may be written in two forms: 


d£ 7P _ yqMo J 

dU U - aj U - aj 


or 


dU 


U 




^ 7 (P - aM 0 2 j) 7 (P - aM 0 2 j) 


W 

(75) 

(76) 

(77) 

(78) 


Equations (77) and ( 78 ) became first-order linear differential equations 
if j is assumed to be a function of U in equation (77) or of P in equa- 
tion ( 78 ). The solution for either equation may be found in most differential 
equation texts (for example, page 116 of ref. 10) . The solutions of equa- 
tions ( 77 ) and ( 78 ) are, respectively, 


f 7 dU 


r 7 an 


p = e J U-af(U) 

r J u-aftu)^ + c 

J U - af (u) 


where j = j Q f(U) and 
. f dP 


r dp 

_ 1 

U = e 7 [P-aMo 2 f(p[| 

/ 

L. 

± f -M J 7[P-aM 0 2 f(P)] dP + C 

7 (p - aMo 2 f(P)J 


where j = j Q f(P). 


(79) 


( 80 ) 
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The constants must he evaluated from the initial or reference conditions. 
The expressions of P as a function of U (eq. (79)) or U as a function of 
P (eq. (80) ) may now he substituted into equation (7*0 or (75) , respectively. 
The resultant equations give U as a function of x if j is a known function 
of U and correspondingly give P as a function of x if j is a known func- 
tion of P. Once either of these variables is determined with respect to x 
and with the use of equations (75a), (75b), and (75c), the density ratio, Mach 
number, and temperature ratio are readily determined. 


Constant Current Density 


The operation of a channel in which the current density and the cross- 
sectional area are held constant is quite practical for experimental use, since 
the driving source required is a constant-current system for each electrode 

used on the channel. If the 



electrodes are evenly spaced 
areawise, the current den- 
sity should be approximately 
constant. 

The solution of the 
equations which represent 
the properties of the flow 
in a constant-area, constant- 
current magnetohydrodynamic 
channel is obtained by sub- 
stituting j = into 

either equation (79) or 
equation (80), integrating, 
and evaluating the constant 
to determine P or U 



(a) Mo = k; 0 < a < 1. 

Figure 25- - Distributions of U, P, M, and T in 
a magnetohydrodynamic channel with B, a, j, 
and A held constant. 


which may then be substi- 
tuted into equation (7*0 or 
(75), respectively. Inte- 
gration of the resultant 
equation then gives U or 
P as a function of x. 

The substitution of 
j = j Q or f (U) = f(P) = 1 

into equations (79) and (80) 
results in 

p - ^ - K 2 - x ) (rtf f 

(81) 


5*1 


and 


U = a - (a - 1) 


P ~ aMp 2 

1 - aMo"2 


-1 h 


( 82 ) 


Equations (8l) and (82) 
may now be substituted into 
equations (74) and (75) , 
respectively. Integration 
of the resultant equations 
and appropriate evaluation 
of the constants give the 
following relationships for 
U and P as functions of 

y: 



and 


y = 


p - i 

7M 0 2 
- (a - !) 


'p - aM 0 2 \ _1//7 
yl -^O 2 / 




Ob) Mo = 0-5; o g a g X. 


(84) 


Figure 23 .- Continued. 


Once U and P are determined as functions of y, then p', T, and M may 
be determined from equations (75a), (75b) , and (75c). Typical values of U, 

P, M, and T are presented in figure 25- The curves are presented for a 
supersonic (M 0 = 4) reference condition and for a subsonic (Mq = 0.5) reference 
condition. 

Several interesting features are seen from the figure. The velocity 
ratio U increases almost linearly for the initial value Mq = 4 provided 
a is less than about 0.9- Also U as a function of y is seen not to be 
sensitive to changes in a; thus, if the conductivity a is sufficiently 
large so that a < 0 . 9 , the velocity will not be sensitive to changes in a. 
However, even though U is not sensitive to cr, it is seen that P, M, and T 
do change a great amount with respect to a and hence are very sensitive to a. 


55 



However, in the subsonic case for a < 1 , it is seen that P and M are 
insensitive to a whereas U and T become sensitive. It may also be 
observed that a range of a exists (fig. 23(c)) for which subsonic velocities 
may be increased, specifically the region 1 < a < 1 /Mq 2. In this region, the 
velocity ratio asymptotically approaches a and does not increase continuously 
as for the supersonic initial condition (fig. 23(a)). 





(c) Mo = 0 . 5 ; 1 < a g 4. 
Figure 23 • - Continued. 


The flow conditions beyond the point of turning represented in figure 23 (a) 
are impossible to attain because the flow cannot reverse its direction as indi- 
cated by the figure. Examination of equation (83) shows that this condition 
occurs for supersonic flow if a is greater than 1 and for subsonic flow if a 
is greater than l/M 0 2 . Since these flows are impossible, the flow will auto- 
matically adjust itself. In the subsonic initial condition such adjustment will 
occur because of upstream pressure changes and readjustment of the entire up- 
stream flow to conform to the condition that M = 1 at the end of the channel. 
If the initial flow is supersonic, a normal shock must occur at some point 
within or ahead of the channel such that the flow will proceed subsonically to 
the end of the channel and emerge at a Mach number of 1 . or less. 
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0 .1 .2 .3 .4 0 .1 .2 .3 .4 

y y 



Figure 23-- Concluded. 


Current Density Proportional to Velocity 

An interesting mode of operation occurs if the current density is assumed 
to vary directly with the velocity ratio U. Since for this mode j = j 0 U, the 
equations are obtained with the substitution of f(U) = U into equations (7*0 
and (79)* After the substitution is made. 



The use of equation ( 85 ) and of j = j Q U in equation (7*0 results after inte- 
gration in 
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Ejp 

Po u o : 


-( x - x o) = y = 


7 + 1 
7 + 1 - a 


log U + 


7 + 1 - a(yM Q 2 + l) 
(7 + 1 - a) 2 ^ 2 


" / 7+l-a \ 

y V 1-a / 


(86) 


Equations (85) and (86) along with equations (73a), (73b), and (73c) may be 
used to determine completely the flow in a magnetohydrodynamic channel in which 
B, cr, and A are held constant and the current density is proportional to the 
velocity. 



(a) Mo = 0.5. 



(b) Mo = 2 . 0 . 

Figure 24 .- Several velocity ratio distributions in a 
magnetohydrodynamic channel with B, a, j/U, and 
A held constant . 


Several typical velocity 
ratio distributions are given 
in figure 24 for a subsonic 
initial Mach number (Mq = 0.5) 
and for a supersonic initial 
Mach number (m q = 2.0). For 
the subsonic Mach n umb er, it is 
seen that if a is less than 
1, the flow decelerates; where- 
as if 1 < a < - 7 - + 1 — (1.77 

yM 0 2 + 1 

where 7 = 1.4), the flow con- 
tinuously accelerates. How- 
ever, although it is not obvi- 
ous from the figure, if a is 

greater than 7 * — — , the 
7M 0 2 + 1 

flow reverses upon itself. 

This flow reversal indicates 
that these initial conditions 
could not occur if the channel 
length is too great. Corre- 
spondingly, for the supersonic 
velocity distributions 
(fig. 24(b)), it is seen that 
if a is less than 1, the 
acceleration is positive and 
the velocity continuously 
increases with y. In this 
case, if a is greater than 1, 
the flow decelerates and 
reverses direction. (See 
curves for a = 2 and a = 3? 
fig. 24(b).) As previously 
discussed, the flow will 
adjust itself to eliminate the 
reversal condition. 


.a.. 
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It is of interest to divide the parametric field of a and M 0 into 
regions showing the values of the parameters that determine acceleration or 
deceleration and direct or reversed flow. This field is mapped in figure 25- 
The nature of each of the indicated regions is determined hy examination of the 
derivative dy/dU evaluated at 


U = 1 and of the conditions 
corresponding to dy/dU = p. The 
first expression 

(a - j \ = 1 - M 0 2 

' dU 'U=l ( a - l)!^ 2 

is used to determine whether U 
increases or decreases at U = 1. 
Examination of the complete expres- 
sion for dy/dU shows that U is 
monotonic with respect to y pro- 
vided flow reversal does not occur; 
thus, if U increases, y will 
always increase or if y decreases, 
U will always decrease for all 
values of y. The term dy/dU = 0 
gives the values of U and y for 
which y is a maximum and hence 
the point at which y reverses 
direction and causes a flow rever- 
sal. This value (u) v of U 

■'max 

is given hy 



(XJ) 


1 max 


1-a 


7 + 1 - a(yM Q 2 + l) 

(1 - a) (7 + 1)M 0 2 


7 + 1 -a 


Figure 25-- Regions showing normal flow and 
reversed flow in a magnetohydrodynamic 
channel with B, a, j/U, and A held 
constant . 


If in a given region U increases and 


decreases and (u). 


(U) v is greater than 1 or if U 
■'max 


is less than 1 a reversal must occur since U can 


'max 


achieve the value (U), 


If U increases and (u)-, 


decreases or if U 


/ w • _L-L U iiiUi CaDCD CX1J.U. YU / -ir 

J'max •'max 

decreases and (U)y mr _ y increases, U cannot achieve (u)y and no flow 


reversal can occur and the flow will proceed through the channel. There are 
two regions in which the term within the brackets of the equation for 

-\r IS 

•'max 

complex and hence is unattainable. The flow in these regions is therefore 
normal. The cross-hatched line indicates the division between the regions in 
which normal flow and reversed flow occur. 


(XJ) 


^max 


is negative, and the exponent is less than 1 so that (U) 
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Constant Electric Field 


A constant-area channel in which the electric field is held constant is 
also quite practical for experimental consideration even though the voltmeter 
reading must he corrected for anode and cathode drops before the field in the 
gas can be determined. However, if the electric field in the gas can be held 
constant, the current density may be expressed either as j = Bcru n (E - U) or 
Jo(E - U) 

j - . The values of U and P as functions of y may be obtained 

E - 1 

with the substitution of f(U) = ^ ^ into equation (79) "to determine P as 

a function of U and substitution of the resultant function into equation (7^) 
to determine U as a function of x or y. When these substitutions are 
made, the results are 

P = y ( y ~ 1 ) m q 2 ( u ~ 1 )( gE - u - i) + g|r - (y - i)e] (8 . 

2 [7U - (7 - De] 


[(7 - l) 2 E 2 - 7F] (U - 1) 

(7 - (7 - i)e] [ 7U - (7 - 1)E] 


L) 2 E 2 - 7F 7U - (7 - l)E (2 - 7)E 2 + F E - IT 

7E ° S 7 _ (7 _ i) E + E ° S E - 1 


where 

_ 7(7 - 1)(2E - l)Mp 2 - 2 [7 - (7 - DE] 

(7 + 1 )Mq 2 

and 

a = (7 - 1)(E - 1) 

and p', M, and T maybe determined from equations (75a) , (75b) , and (75c). 

Plots of P, U, M, and T as determined from equations (87), (88), and 
(75b) and (75c) are presented in figures 26 and 27- A set of values of U, P, 
M, and T are presented in figure 26 for a subsonic (Mo = 0.5) initial condi- 
tion. Examination of figure 26 shows several interesting points. One point is 

that if E is less than 5-5 or y ^ ± ve d°city i s decreased but the pres- 

sure ratios and temperature ratios are increased. Contrastingly, if E is 

greater than "the velocity will approach rather than E 

as might be expected. This point is demonstrated upon examination of equa- 
tion (88) in which it is seen that two terms occur that mate y infinite as 


ko 



U approaches 


(7 - 1)E 


Hence the acceleration attainable from a subsonic 


stream is limited not to E but to 


(7 - 1)E 


This limitation of acceleration 


causes an even more rapid increase in both pressure and temperature. The Mach 
number is seen to decrease for all values o'f E and to be rather insensitive 
to the value of E. 



(c) Mach number. (d) Temperature ratio. 

Figure 26. - Distributions of U, V, M, and T in 
a magnetohydrodynamic channel with B, cr, and 
A held constant, j = BaU 0 (E - U) ; Mq = 0 . 5 . 


Values of U, P, M, and T computed for a supersonic initial Mach 
number (Mq = 2.0) are presented in figures 27 (a) and 27 (b). If the initial 
flow is supersonic, the flow properties are greatly different depending on 


whether E is less than or greater than 


7 

7-1 


or 5.5 when 7 = 1.4. 


kl 



presentation of the results is therefore divided into two sections: (l) E 

with values less than 3*5 (fig* 27(a) ), and (2) E with values greater than 
3*5 (fig* 27 (h)). 




(a) j = BcrU 0 (E - U), E < 3-5, and Mq = 2.0. 

Figure 27 .- Distributions of U, P, M, and T in a 
magnetohydrodynamic channel with B, a, and A 
held constant. 


Examination of figure 27(a) shows that if E is less than — Z — or 3*5, 

7-1 

the velocity increases asymptotically to U = E. Pressures, temperatures, and 
Mach numbers also increase if E is greater than 1. However, if E is less 
than 1 (see curves for E = 0-5), the flow exhibits several irregularities. It 
is seen that the velocity increases indefinitely. Since the channel acts as a 
generator if E is less than 1, such a velocity increase can only be made with 
the reduction of pressure and temperature to negative values. Such values of 
both pressure and temperature are impossible; hence, there must exist some 
value of pressure below which the continuous equations used for this analysis 
no longer apply. It is also interesting to note that the velocity curve for 
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E = 5.^99 indicates a similar mode of transition "between the change from an 

asymptotic approach of U to E for E < — Z — to the approach of U to 1 

7-1 

at E = — 2 — as was observed in figure 21. 

7-1 






Cb) j = BcU 0 (E - U), E > and- Mo = 2.0. 

Figure 27-- Concluded. 


Examination of the velocity of the flow for E greater than 5-5 
(fig. 27 (b)) shows that these flows are forbidden inasmuch as the length 
coordinate y reverses when U attains a certain value. The Mach number 
curves show that this reversal occurs, as in previous cases, when M = 1. It 
is expected that a shock will occur and reduce the velocity to a subsonic value, 
from which it may proceed, as shown in figure 26 , through the channel. 



Constant Velocity 


The constant-velocity condition is treated by letting aj = U; then dU 
of equation (72a) is equal to zero and U becomes constant and equal to 1. 

Hence, aj = aj Q = 1 and since a is constant and equal to j must 

Bcuq 

equal j Q and thus j is also constant. The pressure equation (eq. (72b)) 
may now be written, with these substitutions and m2 = Mo 2 /p as 

dP _ ? M o 2 B jp M 2 - 1 
P dx m 2 - 1 p 0 u 0 2 P 


dP = ?M 0 2 dy ( 89 ) 



(a) Velocity ratio. (b) Pressure function. 



(c) Mach number. (d) Temperature function. 

Figure 28. - Distributions of £— m, and S - . 1 

7Mq 2 7Mq 2 

in a magnetohydrodynamic channel with B, a, A, 
U, and J held constant and BdTJ 0 = (7 - l)j Q . 
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Then, integrating and evaluating the constants gives 

P - 1 + 7M 0 2 y (90) 

From equations (75b) and (73c) 

M = — = M 9 - (91) 

ft \/l + 7Mo 2 y 

and 

T = UP = P = 1 + 7M 0 2 y (92) 


Plots of — - — , T ~ , and M are presented in figure 28 to give a rapid 

7Mq 2 7Mq 2 

estimate of the variations in these quantities. The notable feature is that 
the Mach numbers reduce to subsonic values and approach each other too closely 
to show on the figure after y is greater than about 4. 


Constant Pressure 


The mode of operation in which both the area and pressure are held con- 
stant has several features of interest applicable to the design of practical 
channels. Since the pressure is constant, it should be possible to design a 
constant-area channel having a minimum of wall contact with the plasma. 


This constant -area, constant-pressure mode is created by letting ajM , = U 

in equation (72b). The substitution of this value of U renders dP equal to 
zero and hence P is a constant and therefore equal to 1. Also from equa- 
M„2U _ . UM 2 


tion (73c), M 2 = 


or P = 1 = 


hence 


Since aj equals a constant, j 


or aM 0 2 = 1. This condition is observed to give constant pressure in fig- 
ure 23(a) . The substitution of aj = — , = Mq 2 U, and j = j Q into equa- 

M 2 

tion (72a) results in the following equation for dU: 


M 2 

must equal 


aj = — = -i- = Constant. 


M c 


Mq £ 


j Q , and hence 


<*j = aJ 0 = a = — o 

M 2 


dU 


BJ ° -^ - dy 

Po u o 


(93) 


Integration and evaluation of the constant gives 

U = 1 + y 


( 9*0 
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From equations (73c), 

M = Mq( 1 + y ) 1 / 2 ( 95 ) 

and from equation ( 73 b), 

T = U = 1 + y (96) 

Plots of the values of U, M, and T for various values of Mq computed from 
these equations are presented in figure 29- Mach number, velocity, and temper- 
ature all increase for this case. It is also interesting to note that the 
current density is constant for a constant -pres sure channel as well as for a 
constant -velocity channel. 



(a) Velocity ratio. (t>) Pressure and cur- 

rent density ratio. 




(c) Mach number. (d) Temperature function. 

Figure 29-- Distribution of U, M, J/J 0 j and T in 
a magnetohydr odynami c channel with B, a, P, A, 
and j held constant and BuU 0 = (7 - l)M 0 2 j o . 



An important point must tie observed, however, in case the initial Mach num- 
ber is subsonic. Although equation (94) shows that for a = — i— the velocity 


ML- 


increases indefinitely, it can be shown with the use of equation ( 85 ) that if 


a is the smallest fraction less than 


M, 


2’ 


the flow will exhibit a forbidden 


singularity and therefore completely change its character if the channel is of 
sufficient length. Since it is highly improbable that a can be held exactly 

equal to — — over the entire length of the channel and during the entire 
M 0 2 

operation period, it may well be expected that the increase in U will be 


li m ited, to 


1 

M 0 2 ' 


Constant Temperature 

This channel is often considered as an example of tailoring the current 
distribution to give a desired condition. (See ref. 1.) The constant- 
temperature condition is interesting because the conductivity remains practi- 
cally constant in such a channel. 

The constant -temperature condition is attained by letting aj = ~ — 

7M2 - 1 

Use of this equation in equation (72c) renders dT equal to zero, and hence T 
is constant and therefore equal to 1. Since T = 1, equation (4c) gives 

(7 - l)U 

M = M 0 U and thus aj = M 2[j2 q * ^ow j = j D wlien U = 1; hence, 


7-1 

aj n = - , from which j = 

7Mq 2 - 1 


. _ (yMp 2 


L ) U Jc 


7M 0 % 2 


and a = 


These rela- 


7M r 


- 1 


tions give the necessary tailoring of j to achieve a constant temperature. 

The current density j is seen to vary as U if 7 Mq 2 U 2 j_ s smaller than 

1 and to vary as l/u if 7M G 2 u2 is muc h larger than 1. Typical values of 
this current density variation are shown in figure 30 (d) . 


If these relations for M, aj, and j are substituted into equation (72a), 


the result after several manipulations is 

2 tt 2 


( yM 0 U - dU = 7M 0 2 (7Mq 2 - l)dy 


(97) 


Integration and evaluation of the constant give 


y = 


7M 0 2 (7Mq 2 - l) 


( 7 2 M Q 1 hj 2 + l)(u 2 - l) 

2U 2 


27 M 0 log U 


(98) 
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and from equation (75b) 

P = £ (99) 

Also 

M = MqU (100) 

Plots of the values of U., P, and M as calculated from equations ( 98 ), 
(99), a nd (l00) are presented in figure 30 * The flow properties are rather 


sensitive to Mach number when Mq = — = 0.84515. Also velocity and Mach num- 

\fr 



(a) Velocity ratio. (b) Pressure ratio. 



y y 


(c) Mach number. (d) Current density 

ratio . 

Figure 30.- Distributions of U, P, M, and j/j D 
in a magnetohydrodynamic channel with B, a, 

A, and T held constant. 
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Constant Mach Number 

The current -density distribution may also be varied such that the Mach 
number of a constant-area channel will remain constant. This condition is 

q f + cn* "h-rp *1 -1 — ? ->n i-p /-v A _ (7 + i)u 1 1 'Vi n n 4-Vv -1 n -T rn 1 W m -5 4 -w-1 nrtii n 


attained by letting aj = 


7M 0 2 + 1 


The use of this value for aj in equa- 


tion (72d) makes dM = 0; therefore, M becomes a constant. Also at U = 1, 

7+1 

j = j.., and hence aj^, = a = ^ , and hence j = j-TJ. The substitution of 

° 7Mo 2 +1 

these relations for aj Q and j along with M = Mq into equation (72a) 
results in the following relation between dU and dy: 


dU = 

U 7M 2 


Mn= ® 4 3 2.5 


Integration and evaluation of the constant gives 

2 


7M r 


U = e 


7Mo 2 +1 


(102) 


The values of P and T may he determined from equations (75b) and (75c) and 
are 


7M 2 


_ ml/ 2 _ p ^ M o 2 +l' 


-y 


P = u = = e 


(105) 


Plots of the values of P, U, T^/ 2 , and j/j are presented in figure 31- 


SUMMARY OF RESULTS 


In an analysis of the equations representing the flow of an inviscid plasma 
through steady crossed electric and magnetic fields, the following phenomena 
were found to occur: 

(1) The velocities were limited for certain modes of operation. These 
limitations were usually the result of the balancing of the body forces and the 
pressure gradients. 

(2) The velocity was found to decrease for a number of operational condi- 
tions. In certain of the variable-area modes, this deceleration was sufficient 
to reduce the velocity to zero. 

(3) Flow reversals were observed in constant-area channels when the Mach 
number became equal to 1. Such a situation is impossible; hence, the flow must 
readjust itself to new values. 

(^) Rapid changes in the velocity may occur with small parametric changes 
in channels operated with parametric values near to those for constant velocity. 
The flow in such channels may therefore be unstable. 

(5) The constant-pressure modes of operation exhibited consistent accelera- 
tions, a minimum of the irregular phenomena being found to occur in the constant- 
density modes of the constant-area modes. 

(6) These observations indicate the care that must be exercised in the 
choice of an operational mode as well as of initial conditions if a desired 
result is to be obtained. 


Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va. , May 27, 196^. 
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